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METHOD AND APPARATUS FOR HIGH DATA RATE 
WIRELESS COMMUNICATIONS OVER WAVEFIELD SPACES 

REFERENCE TO RELATED APPLICATIONS 
This application claims priority to Provisional Application No. 60/136,699, 
5 filed on May 28, 1999, incorporated herein by reference. 

FIELD 

The invention is directed to a method and apparatus for high data rate wireless 
communications over wavefield spaces. In particular, the invention provides a high 
data rate wireless communication between transceivers using a multi-dimensional 
10 technique. 

BACKGROUND 

The wireless communication industry is responding to consumer demand for 
high speed data services. The increased Radio Frequency (RF) spectrum acquisition 
costs associated with the achievement of throughputs that are comparable with wired 

15 communication technology has dramatically enhanced the interest in new methods to 
obtain spectrally efficient modulation formats. The technical challenge posed by these 
emerging requirements is residing in the intrinsic physical limitations of the wireless 
propagation environment. Noise, interference, severe multipath and time-varying 
characteristics of the fading channel contribute to make high-speed wireless data 

20 communications a difficult and challenging problem. Traditional communication 
engineering theory states that an increase in data rate is achieved by a combination of 
an increase in bandwidth occupation, an increase in radiated power at the transmitter, 
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and an increase in dimensionality in the signal space (that is the order of the digital 
modulation). A greater description of the prior art and background is provided in 
Appendix 1 of Provisional Application No. 60/136,699, filed on May 28, 1999 
incorporated herein by reference. 

5 Conventional communication technology is depicted in Figures 1 A-B. A 

communication network between transceiver A and transceiver B employs a two- 
dimensional communication protocol using time and frequency. The message to be , 
communicated is broken into a plurality of smaller messages, or packets, that can be 
sent in the various time slots and frequency slots, according to the communication 

10 protocol. Since the receiver knows the communication protocol, it can reassemble the 
packets to replicate the original message. This type of communication protocol is well 
known in the art. A greater description of the prior art and background is provided in 
Appendix 1 of Provisional Application No. 60/136,699, filed on May 28, 1999 

incorporated herein by reference. 
15 what is needed is a communication system that provides high data rate 

communication in a wireless environment. The present invention analyzes practical 
methods to demodulate high data rate signals transmitted from different positions 
(antenna elements), that are afflicted by arbitrarily time-varying fading characteristics, 
with arbitrary time dispersion (frequency selective fading). This is of significant 

20 practical interest because in practice the fading environment can never be considered 
perfectly static even if transmitter and receiver are not in relative motion. The 
presented method and apparatus do not need training sequences or signals to estimate 
the channel: this gives additional advantage in terms of throughput with respect to 
prior art. The method and apparatus is implemented using available hardware for 

25 digital radio transceivers. 

SUMMARY 

The invention overcomes the conventional limitations and provides a high data 
rate wireless communication system between transceivers using a multi-dimensional 
• technique. An exemplary embodiment of the invention for communicating 
30 information from one location to another through a wavefield space using a multi- 
dimensional communication protocol comprises a transmitter configured to receive a 
message and to encode a plurality of signals associated with the message. A first 
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plurality of antennas is coupled to the transmitter and configured to transmit the 
plurality of signals associated with the message into the wavefield space. A second 
plurality of antennas is configured to receive the plurality of signals associated with 
the message. A receiver is coupled to the second plurality of antennas and configured 
5 to decode the plurality of signals associated with the message, wherein the receiver has 
no advance knowledge of at least one of the dimensions of the multi-dimensional 
protocol. The receiver reconstructs the message to retrieve the original information. 

In one aspect of this embodiment, the multi-dimensional protocol is at least a 
three-dimensional protocol and the transmitter includes a plurality of transmit modules 

10 that each encode at least a portion of the message. The receiver includes a plurality of 
receiver modules that each decode complementary portions of the message to 
reconstruct the original information. 

In another embodiment, the receiver is configured to execute a statistical 
procedure to identify and decode the plurality of signals. In one aspect of this 

1 5 embodiment, the receiver is configured to execute a Generalized Likelihood Detection 
(GLD) procedure to identify the plurality of signals and then execute a Maximum 
Likelihood Detection procedure (MLD) to decode the plurality of signals. In another 
aspect of this embodiment, the receiver is configured to estimate the wavefield space 
using a virtual wavefield space that approximates the wavefield space. 

20 In yet another embodiment of a digital wireless transceiver, the transceiver 

comprises a digital circuit for detection of U number of different signals after radio 
frequency intermediate frequency downcon version comprised of digital baseband 
downconversion, trellis search processor, vector symbol synchronizer, wavefield 
estimator, wavefield gradient-based estimator, deinterlever, channel decoder and data 

25 extractor; and a digital circuit for detection of U number of different signals after radio 
frequency intermediate frequency downconversion able to cope with and exploit 
frequency selective multipath propagation, arbitrarily time- varying fading 
characteristics, non Line-Of-Sight communication. 

Advantages of the invention include the ability to communicate high data rates 

30 from one location to another, where the data rates are in excess of those conventionally 
achievable. 
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BRIEF DESCRIPTION OF THE FIGURES 
Figures 1 A-B depict conventional communication technology according to the 

prior art; t 

Figures 2A-B depict communication technology and a three-dimensional 

5 communication protocol according to an embodiment of Reinvention; 

Figure 3 is a flowchart depicting a sequence of events characterizing a 
communication protocol between two wavefield transceivers according to an 
embodiment of the invention; 

Figure 4 illustrates the preferred embodiment of the invention in a transceiver 
(that is an apparatus for transmission and reception of digital information over the 
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wavefield space); 

Figure 5 illustrates in more detail the Radio Frequency (RF) section of the 
transceiver 

with Analog to Digital and Digital to Analog conversions; 

Figure 6 illustrates in more detail the Digital Signal Processing section of the 

Receiver; . 
Figure 7 illustrates in more detail the Digital Signal Processing section of the 

Transmitter; and 

Figure 8 is a flowchart procedure to find the optimum P and M. 

DETAILED DESCRIPTION 
The detailed description of the invention is organized in subsections where 
each subsection explains a topic that is fundamental to explain and justify the methods 
and associated apparatus. This description incorporates Appendix 1 of Provisional 
Application No. 60/136,699, filed on May 28, 1999 herein by reference. 

25 A. Introduction 

The present invention provides high data rate communication between a 
plurality of base stations where each base station transmits, using a plurality of 
antennas, into a multi-path wavefield space. The signals are received by selected base 
station and resolved and decoded to recover the original transmitted informaUon. 

The wavefield space is the space spanned by the channel parameters that 
characterize the multipath fading environment. The communication involves sampling 
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the wavefield space. Sampling the wavefield space creates a multi-link connection 
resource between the transmitter and the receiver. One aspect of the invention 
involves the spatial sampling of the wavefield space by antennas. As the transmitter 
symbols are modulated and simultaneously transmitted using signals that occupy the 
5 same frequency portion of the spectrum, they are distinguishable because of their 
different positions in the wavefield space. The shape of the wavefield space depends 
on the particular propagation environment. The optimum demodulator is a device that 
estimates the wavefield space parameters and performs optimum separation of the 
different signals to obtain the multiple streams. Evidently, the throughput of such a 

1 0 communication system is an integer multiple of the throughput that is achievable at 
each different position in the wavefield space. The concept of wavefield modulation 
can be better understood with reference to Figures 1 and 2, showing the transfer of 
digital information units from the baseband digital signal processing (DSP) domain to 
the Radio Frequency domain and vice versa. Figure 1 depicts the conceptual operation 

15 performed by a traditional transceiver (that is a transmitter-receiver apparatus based on 
prior art). The baseband DSP processors 507 generate digital information units 
(packets) 503 in a non-synchronous fashion. The traditional transceiver 506 is 
responsible for modulating such digital information units (packets) 503 and allocate 
504 each unit (packet) in a preassigned position in the time-frequency two- 

20 dimensional space: one point in the time-frequency plane corresponds to a signal that 
exist in a particular time instant given by the projection of the point onto the time axis 
501 and with an instantaneous frequency given by the projection of the point onto the 
frequency axis 502. The traditional transceiver 506 is also responsible for 
demodulating such digital information units (packets) 503, extract 505 each unit 

25 (packet) from a preassigned position in the time-frequency two-dimensional space, and 
route the units 503 to the baseband DSP processors 507. The traditional transceivers 
are limited in throughput by the number of available slots in the time-frequency two- 
dimensional space since no two signals can occupy the same position in the space. 

In contrast, the wavefield transceiver 506 shown in Figure 6 is capable of much 

30 higher throughput given the same number of slots in the time-frequency two- 
dimensional space since, in fact, two signals can occupy the same position in the 
space. The basic mechanism that allows this is the fact that the wavefield transceiver 
506 is capable of enabling the allocation of digital information units {packets) 503 in a 

5 
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third dimension 601, which is defined as the wavefield space. Basically, the 
throughput is now limited by a three-dimensional volume instead of the traditional 
two-dimensional volume. The wavefield transceiver 506 enables the allocation 504 
(or extraction 505) of digital information units (packets) 503 in the wavefield space 
5 by transmitting (and receiving) from different points in space that exploits the 
multipath propagation environment for communicating across new dimensions. 
Figure 3 describes the fundamental sequence of events 701, 702, 703 that 
characterize a communication session between two wavefield transceivers. 

Given a perfect estimate of the wavefield space parameters, the Maximum 
10 Likelihood principle is the optimum strategy for detection. When the wavefield space 
parameters are not known in advance, the Maximum Likelihood detector can not be 
immediately implemented. A typical approach is the Generalized Likelihood Detector 
(GLD), which uses a Generalized Likelihood statistic (GLS) to detect the signals. The 
Generalized Likelihood statistic, derived for time-varying wavefield spaces, involves 
15 finding the orthogonal projection of the received signal onto the signal subspace; the 
transmitted signal (corresponding to one of the transmitted sequences of symbols) is 
known to lie in a subspace but its exact location is not known because the wavefield 
space parameters are unknown. Detectors in this class are also used in signal 
processing because their statistic is "matched" to the a-priori known signal subspace S. 
20 The (perfectly matched) orthogonal projection of the received signal onto the signal 
subspace S is difficult (if not impossible) to compute. To estimate the subspace S, the 
invention employs virtual subspaces, say V, that approximate the original signal 
subspace S in some sense, and whose orthogonal projections are more easily 
computed. These reduced-size signal subspaces V are obtained by a decomposition of 
25 the fading channel time variations using orthonormal wavelet bases. 

The merits of the wavelet decomposition are to be reconducted to an 
outstanding localization capability in the time-frequency plane. Using expansions for 
a generic signalby means of orthonormal functions that are well-localized in the time- 
frequency plane implies that only a few coefficients of the expansion can be adopted to 
30 represent with extreme accuracy the original signal. Basically, the method and 
resulting apparatus use an efficient representation of the fading statistics in the 
maximum likelihood formulation of the detection problem. 
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In fact, the set of linear vector spaqes defined by such decomposition provides 
also a nested sequence of subspaces that* at increasing level of detail, are "efficient" 
representations of the original signal subspace S. Part of the invention is a method to 
"focus" these (subspace) representations to the original signal subspace using the 
5 concept of Kolmogorov n width for signal classes. The use of the focused wavelet- 
based model in the statistic of the detector results in an effective approach to wavefield 
space demodulation. 

B. Exemplary Transceiver 

The basis for the present invention is communication between a plurality of 

10 base stations where each base station transmits, using a plurality of antennas, into a 
wavefield space. The signals are received by a selected base station and resolved and 
decoded to recover the transmitted information. The invention is performed by a base 
station structure similar to that described in U.S. Pat. No. 5,678,268, assigned to the 
same assignee as this application and incorporated herein by reference. 

15 Wideband Radios use a single Radio Frequency path to access the entire 

allocated band, performing traffic channelization via digital signal processing. 
Traditional radios perform traffic channelization in hardware and typically require an 
additional hardware combining network for multi-channel systems. The flexibility 
increase inherent to wideband stems from the fact that the channelization process can 

20 be done in software or reconfigurable logic which can not only be upgraded but could 
be dynamically changed on a connection by connection basis as described in U.S. Pat. 
No. 5,678,268. Figure 4 shows a Wideband architecture for a transceiver (that is an 
apparatus able to transmit and receive exploiting the wavefield modulation previously 
described). A plurality of receiving antennas 100 are connected to the site RF 

25 distribution 111 which is a network of RF signal distribution circuits. The Power 
Amplifiers 101 are necessary to amplify the transmit signal to reach the remote 
receiver at adequate distance. The Wideband Receiver Modules 102 (one per antenna) 
are in charge of filtering, amplifying, downconverting and digitizing the RF spectrum 
of interest. The digital samples at the output of the ADC from each of the Wideband 

30 Receiver Modules 102 are routed appropriately time-multiplexed to the Wideband 
Digitized Spectrum Data Distribution 108. Digital Receiver Processing modules 105 

7 
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perform traffic channelization, demodulation and data extraction. On the transmit path 
the digital transmitter 106, the Wideband Exciter 103 perform the reverse operations. 

The Wideband Receiver 102 and Exciter 103 RF Modules functional diagrams 
are illustrated in Figure 5. The wideband receiver 102 is designed to be driven from a 
5 site RF distribution 1 1 1 network and minimizes this interface requiring a single, low 
gain RF feed per antenna element. High dynamic range amplifiers 203 and mixers 206 
are utilized to obtain a large instantaneous dynamic range, thus preserving signal 
fidelity. The signal as collected by the antenna 1 00 is filtered 202, amplified 203, 
further filtered 203 and amplified 204. It is then mixed 206 by a first Local Oscillator 
10 209 and filtered 207, amplified 203 and mixed down 206 to intermediate frequency by 
a second Local Oscillator 209. A filter 208 removes unwanted high frequency 
components and the resulting signal is further amplified 203 before being injected to 
the analog to digital converter 210. The clock to the analog to digital converter 210 is 
generated by a frequency source 209 locked to the site frequency reference 109. The 
15 samples at the output of the Analog to Digital Converter are routed to the 108 
Wideband Data Distribution Bus. This Bus 108 is a high speed Time Division 
Multiplexed resource that connects any RF module to any Digital receiver module in a 
non-blocking fashion. All frequency and clock conversion oscillators are locked to the 
external site frequency reference 109 ensuring coherent operation. The Wideband 
20 Exciter 103 RF Module is also illustrated in Figure 5. The samples from the 108 
Wideband Data Distribution Bus are routed to the Digital to Analog Converter 212. 
The analog signal then is filtered 208, amplified 203 and mixed 206 using the second 
Local Oscillator 209 frequency. The RF signal is further filtered 207, amplified 203 
and upconverted 206 using the first Local Oscillator 209 frequency. After proper 
25 amplification 203 and filtering 204 the signal goes to a multicarrier power amplifier 
that is responsible for providing enough RF power to the signal to reach the remote 
receiver after antenna radiation. 

A block diagram for the Digital Receiver Module is shown in Figure 6. The 
traffic channelization (that is baseband conversion) is performed using digital 
quadrature downconversion implemented by means of specialized high-speed ASIC 
301. The samples from the Wideband Digitized Data Distribution 108 are 
downconverted to baseband by means of a quadrature mixer 302 where the mixing 
frequency is generated by a NCO (Numerically Controlled Oscillator) 303. An NCO 
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is a digital implementation of what is known in the analog domain as Voltage 
Controlled Oscillator. The two branches (In-Phase and In-Quadrature) of the samples 
are filtered by filters 314 to remove high freqtiency components of the mixing process. 
Each of the digital downconverters is connected to pulse shaping filters 306 that are 
5 matched to the transmitted pulse shapes and the output of these filters are sampled at a 
rate that is an integer multiple of the symbol rate. These samples are known to 
constitute a sufficient statistic for the trellis search 308 previously described. The 
synchronizer 307 is responsible for optimizing the sampling instants at the output of 
the shaping filter. Wavefield estimation 309 is performed iteratively using the 
10 Generalized Likelihood method described and the gradient update technique 310. The 
processing function following the demodulation 305 are related to the channel coding 
section of the communication system and are well-known to those skilled in the art of 
digital communications. A deinterleaver with channel decoder function 311 follows 
the demodulator. Data extraction 312 refers to the organization of the information 
1 5 bitstream in a way that is compatible with the particular application of the 
communication system. 

The block diagram for the Digital Transmitter Module is shown in Figure 7. 
The data is first properly organized in a bitstream from the user packet data 
distribution bus 1 04 by means of the data extraction function 409. Then, the data is 
20 properly encoded 408 and interleaved 407. The modulator 405 basically organizes the 
high-speed single bitstream in multiple lower speed bitstreams. These bitstreams are 
modulated according to the particular digital modulation format of interest (for 
example M-QAM or M-PSK) and routed to the pulse shaping filters 406. The device 
that follows is a digital up-converter 401 that translates to a first intermediate 
25 frequency the baseband signal at the output of the pulse shaping filters. The digital 
upconverter is constituted by filters 404, mixers (multipliers) 402, a Numerically 
Controlled Oscillator 403 and an adder 41 1. The digital samples at the output of the 
digital upconverter are routed by means of the Wideband Digitized Spectrum Data 
Distribution Bus 108 to the respective Wideband Exciter RF module 103 where they 
30 are converted to analog signals and upconverted to RF. 

Each module 102, 103, 105, 106 is controlled by a respective microcontroller 
215, 216, 313, 410, each of which also provides an interface with the data distribution 
resources (Wideband Digitized Spectrum Data Distribution 108 and User Data 
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Distribution 104). Multiple modules can be replicated to meet traffic capacity and/or 
multiple antenna requirements. 

C. Detailed Description of the Signal Processing Problem 

Consider the communication between two base-stations where the transmit 
5 base-station has U antennas and the receive base-station is equipped with K antennas. 
At each antenna a different symbol stream is transmitted, so that the rate at which 
' symbols are transmitted by the U-antenna transmitter is U times the rate at any single 
antenna. 

It will be immediately evident to those skilled in the art of wireless receiver 
10 technology that the effect of a time-varying multipath channel can be represented at 
baseband (at the output of the pulse shaping filters 306 in Figure 6) using the complex 
sequence of samples 

y k (n) = a(n) T h k (n) + Ti k (n), n=0,l,2,...RN-l, k=l,2 K (1) 

where T means transpose for vectors or matrices, 
15 a(n)=[a,(n),... a,(n-D), a^n),... a^n-D), .... a^n),... a y (n-D) ] T 

is a vector that contains the N symbols transmitted by each of the U antennas 
interlaced (or zero-padded) with R-l zeros and 

h k (n)=[g k>1 (n,0), g M (n,l) gk.,(n,D), g w (n,0), 

g u (n,l) gt2(n,D) g k .u(n,0),...,g k . u (n,D)] T . 

20 is a vector of complex samples that contains the baseband discrete-time impulse 

responses from any of the U transmit antennas to the kth receive antenna. The generic 
element gyjn,m) of h k (n) represents the response of 

1. pulse shaping at the 1th transmit antenna; 

2. transmit RF filters at the 1th transmit antenna; 

25 3. mth tap of the discrete time multipath impulse response at time n between 1th 
transmit antenna and kth receive antenna; 
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4. receive RF filters at the kth receive antenna; and 

5. pulse shaping at the I-th transmit antenna. • 

D is an integer that characterizes the length of the finite impulse response (FIR) 
g w (n,m) between any of the U transmit antennas to any of the K receive antennas. The 
5 index n simply describes the fact that such impulse responses are time-varying. The 
complex scalar ri k (n) represents additive White Gaussian Noise afflicting the 

observation at the kth receive antenna. 

More details regarding this technique of modeling multipath can be found in 
M. Martone, "An adaptive algorithm for adaptive antenna array low-rank processing in 

10 cellular TDMA base-stations," IEEE Trans. Commun., vol. 46, pp. 627-643, May 
1998; M. Martone, "On MMSE Real-Time Antenna Array Processing Using Fourth- 
Order Statistics in the US Cellular TDMA System," IEEE Journal Sel. Areas in 
Comm., vol. 1 6, N. 8, pp. 1 396- 1410, Oct. 98; and M. Martone, "Cumulant-based 
adaptive multi-channel filtering for wireless communication systems with multipath 

15 RF propagation using antenna arrays" IEEE Trans. Veh. Tech., Vol. 47, no. 2, pp. 377- 
391 May 98. 

The samples y k (n) are samples of the baseband signal collected at the kth 
antenna at sampling rate equal to R/T if we denote 1/T as the rate at which the U „ 
antennas transmit synchronous symbols. More explicitly 
20 1 . aj(n) is a complex number representing a point (symbol) in the constellation if n is 
equal to 0, R,2R,3R,...,NR; and 

2. aj(n) is equal to zero if n is not equal to 0, R,2R,3R,...,NR. 

The integer parameter R is the so called oversampling factor and it represents the 
number of signal samples that are contained in one symbol period T. It is well-known 

25 that R=l is sufficient for optimum detection. However in a practical receiver R>1 
makes implementation more robust with respect to synchronization issues. It is 
important to remind that in the model at hand for K=l, L=l, R=l and static channel, 
that is g kJ (n,m)= g,.,(m)= g(m) for any n, equation (1) collapses into the well-known 
model for reception of digital signals in presence of Intersymbol Interference (see J. G. 

30 Proakis, "Digital communications," McGraw-Hill, Third Edition, 1995). 

y,(n) = S m g(m) a,(ii-m)+ Ti,(n), n=0,l,2,. . .RN-1 (2) 

11 
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Equation (2) represents an unrealistic model in real life wireless communication 
channels and much more importantly the communication system modeled by equation 
(2) has a maximum throughput of N symbols per NT seconds, while equation (1) is 
representative of a communication system where the throughput is UN symbols per 

5 NT seconds. 

The task of the receiver base-station signal processing is to obtain an estimate of 
the vector of symbols a(n) for n=0,l, . ..NR-l, (in other words an estimate of the 

sequence a=[a(0) T , a(l) T a(NR-l) T ] T ), g*en only the received samples y k (n). 

If the communication system uses Q-ary symbols then the vector a can be any of 
10 Q m different vectors. We identify the hypothetical vector as s=[s(0) T , s(l) T s(NR- 

1)T. 

The Maximum Likelihood formulation of the detection/estimation problem can be 
stated as in the following: 

15 Find the minimum of 

Ms] = S^....* w^-. iyx< n >- S ^ T h '< n >i 2 

over "any possible" 

s = [s(0) T s(l) T ,...,s(NR-D T ] T 

Observe that "any possible" refers in the above formulation to the fact that since 
20 the elements of s are symbols that can assume only a finite number of values (the 
Alphabet of the digital modulation method of choice), there will also be only a finite 

number of vectors s. Also, one of the many s= [s(0) T ,*l) T s(NR-l) T f_is indeed 

the transmitted one vector a= [a(0) T a(l) T ..,a(NR-l) T ] T . 

However the big practical problem in this optimum formulation of the problem is 
25 that the receiver does not know the multipath channel parameters (the wavefield 
parameters). Prior art has solved this problem sending training sequences of known 
symbols that allow channel estimation. This technique has two drawbacks: 

1. it works only for static channels (that is non-time varying); and 

2. it dramatically reduces the throughput of the system because of the periodical 
30 training. 

The first modification is to formulate the problem differently as in the following: 
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Find the minimum of 

I^fs] « Swi^S^aw |y k (n)-s(n) T h k (n)| 2 

over "any possible" 

s=[s(0) T s(l) T ,...,s(NR-l) T ] T 

5 and over h^n). 

This problem is solved in two stages: 

Step 1 : For any possible sJFind the Maximum Likelihood Estimate of h k (n), call it 
hk[s,n], 

Step 2: Find the minimum of L^s] = 2j kHi2 K_2n=<u,...NR.i lyk(n)-s(n) T h^n]) 2 

1 0 over any possible s. 

This two-step particular solution of the problem is called Generalized Likelihood 
solution (see L. Van Trees/'Detection, Estimation, and Modulation Theory" Part I, 
John Wiley 1968). Unfortunately this way of solving the detection problem still 
places a formidable computational challenge, essentially because a large number of 

1 5 free parameters are contained in h k (n), for any n and any k. In fact h k (n) is a non- 
parsimonious (in other words it has significant redundancy) representation of the 
wavefield multipath propagation environment. To minimize the amount of 
redundancy in h k (n) and get to the possible minimum number of parameters that truly 
characterize h k (n) especially as connected to our detection problem we use a wavelet 

20 orthonormal basis expansion of each element of h k (n), the wavefield response. 

D. Multiresolution Decomposition of the Wavefield Multipath Response 

An unknown deterministic or random signal can be decomposed according to its 
resolution (components at different levels). An excellent and very intuitive 
explanation of the concept of decomposing a signal at different resolutions is 
25 contained in the book "Wavelets and Filter Banks" by G. Strang and T. Nguyen, 
Wellesley-Cambridge Press, 1996. 

The crucial point of almost all of the proposed methods for the detection of a 
signal of the form in (1) (even in the case of a single transmit antenna) is based on the 
assumption that the multipath channels are quasi-static, that is time invariant over the 
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length of the transmitted frame. Slow variations of the channels are then compensated 
by using adaptive algorithms that ultimately force the estimate to be constantly in 
search of a convergence point. If the channel coefficients variations in time are fast 
with respect to the convergence time of the adaptive algorithm, significant degradation 

5 may result. A more reasonable alternative is tp model the time-varying components of 
the impulse response of the multipath channel as low-pass Gaussian processes with 
their correlation representing ensemble average characteristics. These average 
characteristics are in turn used to fit parametric models so that a bank of Kalman filters 
can be used for tracking. This approach is impractical because the parametric models 

10 that best fits the channel dynamics have a large number of free parameters. Moreover, 
the generalization of this approach to a multi-antenna environment where possibly 
multiple different channel responses have to be estimated/tracked, is definitely a 
formidable task in terms of computational effort. 

This invention departs from these well-known and obvious approaches in search of 

15 new methods that more "economically" can characterize the time varying nature of the 
detection/estimation problem. A multiresolution representation of the fading process 
in each component of the channel response is introduced, elaborating some ideas of (J. 
Zhang and G. Walter, "A wavelet-based KL-like expansion for wide-sense stationary 
random processes" IEEE Transactions on Signal Processing, vol. 42, no. 7, pp. 1737- 

20 1745, July 1994); (R. W. Dijkerman and R.R. Mazumdar, "Wavelet representations of 
stochastic processes and multiresolution stochastic models" Transactions on Signal 
Processing, vol. 42, no. 7, pp. 1640-1652, July 1994), M. I. Doroslovacki and H. Fan, 
"Wavelet-based linear system modeling and adaptive filtering," Transactions on 
Signal Processing, vol. 44, no. 5, pp. 1 156-1 167, May 1996. The unknown channel 

25 time variations are in fact decomposed using optimal unconditional bases such as 
orthonormal wavelet bases. 

It is extremely important to realize that modeling of linear systems by basis 
functions can turn a time-varying system identification problem into a time-invariant 
one. Essentially inspired by Fourier harmonic analysis ideas, simple exponential bases 

30 were used in G. B. Giannakis and C. Tepedelenlioglu, •'Basis-expansion models and 
diversity techniques for blind identification and equalization of time-varying 
channels," Proceedings of the IEEE, vol. 86, no. 10, pp. 1969-1986, Oct. 1998. 
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The exponential basis function of a Fourier-like decomposition has infinite 
duration, so clearly any representation of a time-localized signal is not efficient and 
adequate. To connect this consideration with a simple and intuitive example, consider 
a signal made of two pure oscillations occurring at non-overlapping distinct time 
5 intervals. 

The Fourier transform reveals the presence of the two tones, but with no 
localization in time, which forces a wrong representation for the signal. In other 
words, representations based on exponentials lack "parsimony." The time-varying 
reflections of electromagnetic waves caused by a moving transmitter (or receiver) are 

1 0 in nature very similar to a signal where multiple propagation modes are present at 
different times. To represent the frequency behavior of a signal locally in time, the 
signal should be analyzed by functions which are localized both in time and frequency, 
for instance, signals that are compactly supported in the time and frequency domains. 
In reality no function can be compactly supported simultaneously in time and 

15 frequency so the scientific community has directed intense research efforts in the 
development of functions with "good" time-frequency localization: wavelets are the 
best known tool for linear time-frequency analysis. An obvious solution to the time- 
localization problem of the Fourier analysis is to localize the tones in the Fourier 
representation by windowing several time-consecutive intervals of the analyzed signal. 

20 This is in fact the short-time Fourier transform which basically operates a partition (or 
'tiling") of the time-frequency plane in rectangles of equal area. Wavelets can offer a 
different and valuable compromise: the frequency localization is logarithmic in 
frequency. Looking at the time-frequency plane, this means that time localization gets 
finer at the highest frequencies. The wavelet transform replaces the Fourier 

25 transform's sinusoidal waves by a family of functions generated by translations and 
dilations of a single window called a wavelet. Complicated signals can be represented 
using only a few wavelet and scaling functions and statistical signal modeling and 
processing methods based on the wavelet transform are in most cases much more 
effective than classical time-domain or frequency-domain approaches. 

30 We give a brief explanation of the main concept disclosed in this invention by using 
an oversimplified continuous time single antenna model. Consider a sequence of digital 
symbols collected in the vector a, digitally modulated to form the signal waveform x(t,a). 

15 
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The signal x(t, a) is the input to h(t,-c) the time-varying impulse response of a multipath 
channel, as 

y(t) = J t »^., h(t,t-t) x(x,a) dr: 
Assume that the kemelh(t,x) can be represented by an expansion of the form 
5 h(t,r) = S l6 „w,(x)T l (t) 

where {¥,(0) l6 „ is a set of functions (wavelets) and {w,(T)} Ie „ are the (wavelet) 
coefficients of the expansion. Using this representation we can express y(t) as 



y(t) 



- Si.h *i(0 ., w,(t-t) x(x,a) dx = S I6 h ^.(t) x*. (x,a) 



where x', (x,a) = J t _. t w,(t-x) x(x,a) dx. This representation can be economical if the 

10 set {¥,(0} Ie „ contains a small number of functions. The multiplication Y,(0 x\ (x,a) 
emphasizes the time selection of the input content and its frequency translations. 
Observe that x\ (x,a) for any I in the set H is the output of a time-invariant linear 
filtering operation: it is simply obtained filtering x(a,t) by w,(t). Another interesting 
property of this representation is that it gives us the ability to select for any given time 

15 interval only those wavelet coefficients responses w,(t) that affect the output during the 
time interval of interest. The use of wavelet bases allows the selection of the subset in 
H which will model the local and global characteristics of the channel dynamics. A 
clear attempt of the engineering community to have a global characterization of the 
time-varying nature of a multipath propagation link is the well-known WSSUS (Wide 

20 Sense Stationary Uncorrected Scattering) model (J. G. Proakis, "Digital 

communications » McGraw-Hill, Third Edition, 1995)). The WSSUS assumption is 
only an approximation of the real propagation mechanism in a dense multipath 
environment. Indeed h(t, x) has smooth (low frequency) variations with occasional 
rapid changes that can not be possibly modeled as stationary. In particular a channel 

25 model which is uniquely based on ensemble average characteristics will inevitably 
loose local characteristics that are caused for example by rapid transitions and sudden 
changes of the impulse response dynamics (for example the mobile transmitter coming 
to a sudden stop), while the multiresolution model will retain such information in a 
few high resolution "detail" coefficients. The idea is that the wavelet-based 

30 representations can better model not only non-stationary but also stationary processes. 
We will however use the WSSUS model to tune our channel representation and to 
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eliminate gross redundancy: in other words second : order statistics will provide general 
guidelines on how to select the low-resolution coefficients. It is however understood 
that retaining only the low-resolution coefficients will not give us the kind of 
"parsimony" that we expect from a wavelet-based representation. To achieve this goal 
5 a few high resolution coefficients empirically selected will be maintained. 

Going back to our original discrete-time problem, we would like to obtain a 

practical basis representation for the discrete-time channel g^^k) that is easily 
parameterized, it retains the essential features of the fading process and it characterizes 
these features using a small number of time-invariant coefficients. These requirements 

10 can be satisfied using a discrete-time wavelet series representation. The orthonormal 
expansion we will derive can be thought more properly as a multiscale sub-band 
decomposition using filterbanks of the sampled response of a multipath channel. We 
direct the interested reader to P. P. Vaidyanathan, "Multirate Systems and Filter 
Banks," Prentice Hall, 1 993, and in particular to O. Rioul, "A discrete-time 

15 multiresolution theory," IEEE Transactions in Signal Processing, vol. 41, no. 8, pp. 
2591-2606, Aug. 1993, where a "discrete-time" approach to multiresolution 
representations is described. To simplify notation consider g(n,k) as the generic 
response for a generic set of indices il,i2. It is fundamental to observe that g(n,k) can 
be represented at lower resolution applying a half band low-pass filter having impulse 

20 response c 0 (n) followed by a down-sampling-by-two operation 

C,^(k)=S j c 0 (j)g(2m-j,k). 

This equation represents a mapping from a vector space to itself. An added 
"detail" can be obtained by using a high-pass filter with impulse response c,(n), then 
again down-sampling by two 
25 ^W-SjCjCj^g^-jJc). 

The filters Cj(n), 1=0,1 satisfy some regularity constraints and form an orthonormal 
set. If c,(n) (in the z-domain Q(z)) for I = 0,1 is a dyadic perfect reconstruction filter 
bank (P. P. Vaidyanathan, "Multirate Systems and Filter Banks," Prentice Hall, 1 993), 
one can think of C Ilin (k) and ? lJB (k) as a decomposition of g(n,k) onto a discretetime 
30 orthonormal basis with the following (perfect) reconstruction rule 
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g(n,k) = S m C,^(k) c 0 (2 m-n) + S ro $ ljn c,(2 m-n), 

which is indeed a sum of orthogonal projections. Using the same filters Cj (n), 1=0,1, 
the decomposition could be performed again in lowpass and highpass down-sampled 
components, and then further the lowpass component split again, up to a desired 
5 "depth" or resolution, so that one can always "explode" the lowpass component 

resolution 

C jJB (k) = S m C j+1/m (k) c 0 (2 m-n) + S. ^ c,(2 m-n) 

For example consider the equation above for j=l to get 

g(n,k) = S m [ S ra . C 2 . m (k) c 0 (2 m'-m) + S m .^ ra . (k) c,(2 m'-m) ] c 0 (2 m-n) 
10 +S m § ljn (k)c,(2m-n) 

= S m C 2jm (k) Co^ (2 2 m-n) + ZU S ra 5„ (k) c»<2 ' m-n) 
where c 0 <2) (n), c, (J) (n) and c, m (n), 1=1,2 are filters obtained in the z-domain as 
C 0 w (z)=C 0 (z)C 0 (2 2 ),C/ ,) (z)=C 1 (2), C/ 2) (z)= C.CZ 1 ) C 0 (z). 

In the last equation we have used the Noble identity (P. P. Vaidyanathan, 
15 "Multirate Systems and Filter Banks," Prentice Hall, 1993) which says that a 

downsampler by two, followed by a filter F(z) is equivalent to the filter F(^) followed 
by the downsampler by two. 

We can think of the derived expansion for g(n,k), as a decomposition of g(n,k) 
onto a discrete time orthonormal basis at resolution level P=2. The method can be 
20 applied recursively to obtain at generic resolution depth P a filter bank with P+l 
branches. 

The channel parameter g ili2 (n,k) at any fixed k can be considered a signal at 
resolution 0 with an associated signal subspace V 0 = S. Of course gi I42 (n,k) for 
il=lj2 Uandi2=l,2,...,Kandn=0,l,...,NR-l spanning V 0 perfectly represents the 
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wavefield response, but it has high degree of redundancy and it complicates our 
detection problem. 

The orthonormal basis c p (q) (2 q m-n) for p=0, 1 and q=l ,2,. . .P can be used to 
decompose g iKi2 (n,k) at resolution P 

5 gn.i2(n,k) - Z^.^w Xp^mmOO c 0 (?) (2 P m-n) 

+ E MXJ> K^i^ui - 5uuu2(k) m-n) (3) 

The smaller P the higher the resolution (is the highest resolution). So greater P 
means less redundancy in the wavefield response representation. The filters c { (n) for 
i=0,l satisfy some regularity constraints, and form an orthonormal set: they are 
1 0 defined as a dyadic perfect reconstruction filter bank ("Wavelets and Filter Banks" by 
G. Strang and T. Nguyen, Wellesley-Cambridge Press, 1996). We are not interested in 
the details of why these constraints are related to the discrete time wavelet transform, 
all we care is that Cp^i^OO and Cuui-aOO are <low redundancy> representations of 
g ilU2 (n,k). Equation (3) can be rewritten in vector form as 

15 ai lB (n > k)*c(P i n) T w lliB (P > k) 

In particular, given a particular P we can truncate the second summation by 

considering that we can neglect the last 2^ , M NR T m components of the vector 

w a aOO (because typically very close to zero in real life multipath channels) and of the 
vector c(P,n) so we can obtain a parsimonious representation of the channel as 

20 h k (n)= C(P,M) w k (P,M) (4) 

where C(P,M) is simply a matrix of vectors c(P,n) T where the last 2 n « u . j*NR 2 m have 
been eliminated and w k (P,M) is simply a long vector obtained by the concatenation of 
vectors w iU2 (P,k) for il=l,2,...,U and i2-l,2,...,K and k=0,l,...,D where the last 

S^, j^NR 2* m have been eliminated also. 
25 Now that we are equipped with the representation (4) of the wavefield response we 
go back to equation (1) and observe that we can write 
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y k (n) = a(n) T C(P^)w k (P,M) + ti k (n) = a(n) T w t (P,M) + Ti k (n), (5) 

where we have indicated a(n) T = a(n) T C(P,M). 

Denote now d(n) T = s(n) T C(PJvl) and d=[ d(0) T , d(l) T , . . ., d(NR-l) T ] T , and re- 
formulate the Generalized 2-step Likelihood problem as: 
5 Step 1 : For any possible d find the Maximum Likelihood Estimate w k (P,M), call it 
w k [P,M,d], and 

Step 2: Find the minimum of L^d] = S MA ....ic_ S n-o.i....NR.i lyk( n ) " d ( n ) T 

w k [P,M,d]| 2 over any possible d. 

The problem is now in simple form because the Maximum Likelihood estimate of 
10 w k (P,M) given s is dramatically simplified by the fact that the wavefield response is 
characterized by a static vector w k (P,M). 

E. Focusing the Generalized Energy Detector 

This subsection describes a method to select the proper integer parameters M and P, 
given a certain accuracy of representation of the wavefield space. 

15 Described above is the practice that a multiresolution decomposition of the wavefield 
response dynamics obtains a nested subspace representation for the signal 

S=V 0 cV,c ...V, ... (6) 

at increasingly high level of detail. Of course one is interested in a decomposition which 
is the least possible computationally expensive, that is the "smallest" representation 

20 subspace, for a "reasonable" approximation accuracy. 

If the unknown channel variations are characterized statistically (for example using 
the well-known WSSUS, Wide Sense Stationary Uncorrelated Scattering assumption) S 
is not anymore a linear subspace but an ensemble of linear subspaces and it is possible to 
design a representation subspace that is close "on average" to S. 

25 A measure of the distance between two closed subspaces S and V is the deviation 
d(S,V) which is based on the concept of (constrained) Kolmogorov n width. 

The deviation is a one-sided measure. That means that it addresses only the question 
of how well S is included in V but not vice versa, which is also an important issue. In 
fact, V may be too "large" to be of any practical interest. What a practical 

20 



WO 00/74264 



PCT/US00/12069 



implementation really needs is the "smallest" subspace that minimizes the deviation. A 
"small" subspace implies that the number of coefficients needed to represent the 
orthogonal projection is small, which results in less degrees of freedom, reduced 
sensitivity to noise and simplicity of computation. Since the subspace representation for 
5 S over a set of predetermined basis vectors indexed by the integer P (and for each P 
approximated by M) is of interest, the procedure is as follows: 

1. find the deviation for each subspace with respect to S in the class of subspaces 
represented by Equation (6); 

2. select the one (that is select the resolution depth P) that yields a tolerably small 
10 deviation; and 

3. "shrink" the subspace retaining only part of the coefficients and select M in a way that 
the deviation is still close to the one selected in step 2. 

If the unknown channel variations are characterized statistically, S is not anymore a 
linear subspace but an ensemble of linear subspaces. Statistical characteristics of the 
1 5 fading channel are typically available from experimental measurements in the form of a 
scattering function or autocorrelation function. This knowledge translates to the 
knowledge of the following correlations 

E{g n , i2 (n,kl)* g i3 ^(n-^,k2) } =R iUW3<i4 (rn,kl,k2), 

for any il,i2,i3,i4,m,kl,k2. 
20 The average deviation between S and V is 

d(S,V) 2 = E i -Trace[A T ' P A A RJ (7) 

where 

1. Trace[ ] is the Trace operator for a matrix, that is the sum of the elements along the 
diagonal; 
25 2. A= I K ® (Inr® s); 

3. P A = B (B T * B)-» B r ' with B=[d(0), d(l) d(NR-l)]; 

4. is the Kronecker product operator between two matrices (the element-wise product); 

5. * means complex conjugation for scalars, vectors and matrices; 
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6. R„= E{ hh T *} with h= [h,(0) T h,(l) T h,(NR-l) T h 2 (0) T ,.., h 2 (NR-l) T , 
hK(0) T ,...,h K (NR-l) T ] T ;and 

7. E{ } means statistical expectation. 

Detailed and rigorous derivations that explain equation (7) are contained in 
5 Appendix 1 of Provisional Application No. 60/136,699, filed on May 28, 1999 
incorporated herein by reference. 

A technique to obtain the best values for P and M is then 
1. Find the deviation for each subspace with respect to S in the class of subspaces 
represented by Equation (6) over the range 1 i: p £ Pma,,; 
10 2. Select the smallest p such that dCS.V^S/vV,) 2 s A; ™ d 

3. "shrink" the subspace retaining UK (D+l) NR (1- S^,^ 2 m ) coefficients in the 

formation of the function L^d] for M^,, s M £ Mp.- 

Observe that A is a preselected small number that determines the 
approximation accuracy of the representation subspace. 
15 The algorithmic procedure to find the optimum P and M is detailed in the flow 

diagram of Figure 8. 

F. Real-Time Implementation 

There is a celebrated method known as the Viterbi algorithm (G. D. Forney, 
"Maximum Likelihood Sequence Estimation of Digital Sequences in the Presence of 

20 Intersymbol Interference," IEEE Trans. Inform. Theory, vol. IT-1 8, pp. 363-378, May 
1972) that performs the task of minimizing L^d] efficiently. The significance of the 
Viterbi algorithm is that the number of computations required for the detection of a 
sequence of length N grows only linearly with N. The Viterbi algorithm is well 
known to those skilled in the art of Maximum Likelihood Sequence Detection (or 

25 Estimation). 

To apply the Viterbi algorithm one has to reduce the function L^s] to a sum 
of terms where each term corresponds to the length of a branch in a suitable trellis (or 
tree) diagram. Also, one needs a real-time algorithm to compute the the estimated 
vector w k [P,M,d] in LJd] on a sample by sample basis. This can be easily derived 
30 defining w tn [P,M,d] as w k [P,M,d] estimated at time step n. Possible approaches are 
gradient-based methods (like the LMS, Least Mean Squares) and Recursive Least 
Squares-based (RLS) methods. 
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To reduce the function LJd] to a sum of terms where each term corresponds to 
the length of a branch in a suitable trellis (or tree) diagram, it is convenient to define 
the variables: 

F = [d(n) T ,d(n-l) T ,.-,d(n-D) T ] T for n=D,D+l,..., NR 
5 and the quantities: 

Gjhj ( F D , F n+ |)= ^k-l-JC ^m-l-d yk( m ) " °-( m ) T W k.mP > »^» ^ I* 

where w^M, d] is recursively computed as: 

„[P,M, d]= w^JP.M.d] + u [ y k (m) - d(m) T w^.fP.M, d] ] d(m) T m=l,2,...,n, 



and the initialization is w^JP^M, d]=0. 
10 Observe that G n+I ( F D , F n+1 ) depends on d(m), for m=l,2,..,n, that is, on a path in 

the trellis. 
Also, defining 

J,., ( F„, F n+I ) = | y k (n) - d (n) T w^tP.M.d] | 2 

we see that J n+1 ( F n , F.+, ) depends on d(n) and d(n-l) only, that is, on a branch in the 
15 trellis. It is easy to verify that: 

( F D . *V.) - G„ ( F D , .... FJ + J„+! ( F n , F n+1 ). 

There is a one-to-one correspondence between the sequence of transmitted 
symbols and the sequence of states F D , .... F^,. Therefore the problem of selecting the 
most likely symbol-sequence is equivalent to that of selecting the most likely sequence 
20 of states. The Maximum Likelihood detection problem has been reduced to the 

selection of a path through the trellis, once the branches joining states F„ and F n+I have 
been assigned the values taken by the function J rt ( F n , F n+l ) usually referred as the 
branch metric. Any particular digital modulation scheme, with its own associated 
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alphabet, length D, and number of transmit antennas U will generate a different trellis 
structure. The minimum-metric path corresponds to the most likely sequence of states 
and hence to the most likely sequence of symbols. A summary of the necessary steps 
is given below. Step 1 Initialize!: Let n=D, compute G D (F D ) for every possible value 

5 of F D and store them. 

Step 2 Extend survivor paths : Let n=n+l and compute J n+I (F„, F^,) for each 
pair of states F„, F^, such that the transition from F 0 to F^, is allowed by the trellis 
structure. Step 3 Drop Paths: For each state F^, compute iv, (F^) as the minimum 
overanyF^of iv, (F n+ |)+ J„+i (F„» F n+ |) where u D (F D )=G D (F D ) by definition and 
10 iv, (F^,) is recursively computed. Step 4 Trace-back procedure. Step 5 Iterate: Go 

to Step 2 if not n=N. 

The detected symbol is given at the output by means of the trace-back 
procedure, the procedure that outputs the estimated symbol as the one associated with 
the minimum euclidean distance path among the retained paths (also defined 
15 survivors). 

G. Conclusion 

The invention provides many advantages over known techniques. Advantages 
of the invention include the ability to communicate high data rates from one location 
to another, where the data rates are in excess of those conventionally achievable. 
20 Having disclosed exemplary embodiments and the best mode, modifications 

and variations may be made to the disclosed embodiments while remaining within the 
scope of the invention as defined by the following claims. 
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CLAIMS 

1 . An apparatus for communicating information from one location to another 
through a wavefield space using a multi-dimensional communication protocol 
comprising: 

5 a transmitter configured to receive a message and to encode a plurality of 

signals associated with the message; 

a first plurality of antennas coupled to the transmitter and configured to 
transmit the plurality of signals associated with the message; 

a second plurality of antennas configured to receive the plurality of signals 
10 associated with the message; and 

a receiver coupled to the second plurality of antennas and configured to decode 
the plurality of signals associated with the message to reconstruct the original 
information, wherein the receiver has no advance knowledge of at least one of the 
dimensions of the multi-dimensional protocol. 

15 2. The apparatus of claim 1 , wherein the protocol is at least a three-dimensional 
protocol and wherein: 

the transmitter includes a plurality of transmit modules that each encode at 

least a portion of the message; and 

the receiver includes a plurality of receiver modules that each decode 
20 complementary portions of the message to reconstruct the original information. 

3. The apparatus of claim 1, wherein: 

the receiver is configured to execute a statistical procedure to identify and 
decode the plurality of signals. 

4. The apparatus of claim 3, wherein: 

25 the receiver is configured to execute a generalized likelihood detection 

procedure to identify the plurality of signals and then execute a maximum likelihood 
detection procedure to decode the plurality of signals. 
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5 The apparatus of claim 4, wherein 

the receiver is configured to estimate the wavefield space using a virtual 
wavefield space that approximates the wavefield space. 

6. An apparatus for communicating information from one location to another 
5 though a wavefield space using a multi-dimensional communication protocol 
comprising: 

a first transceiver configured to receive a message and to encode a plurality of 
signals associated with the message; 

a first plurality of antennas coupled to the first transceiver and configured to 
10 transmit the plurality ofsignals associated with the message; 

a second plurality of antennas configured to receive the plurality of s.gnals 
associated with the message; and 

a second transceiver coupled to the second plurality of antennas and configured 
to decode the plurality ofsignals associated with the message to reconstruct the 
15 original information, wherein the second transceiver has no advance knowledge of at 
least one of the dimensions of the multi-dimensional protocol. 

7. The apparatus of claim 6, wherein the protocol is at least a three-dimensional 

protocol and wherein: 

the first transceiver includes a plurality of transmit modules that each encode at 

20 least a portion of the message; and 

the second transceiver includes a plurality of receiver modules that each 
decode complementary portions of the message to reconstruct the original information. 

8 . The apparatus of claim 6, wherein: 

the second transceiver is configured to execute a statistical procedure to 
25 identify and decode the plurality ofsignals. 

9. The apparatus of claim 8, wherein: 

the second transceiver is configured to execute a generalized likelihood 
detection procedure to identify the plurality ofsignals and then execute a maximum 
likelihood detection procedure to decode the-plurality of signals. 
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10. The apparatus of claim 9, wherein 

the second transceiver is configured to estimate the wavefield space using a 
virtual wavefield space that approximates the wavefield space. 

11. A digital Wireless Transceiver comprising: 

5 a digital circuit for detection of the U number of different signals after radio 

frequency intermediate frequency downconversion comprised of digital baseband 
downconversion, trellis search processor, vector symbol synchronizer, wavefield 
estimator, wavefield gradient-based estimator, deinterlever, channel decoder and data 
extractor; and 

10 a digital circuit for detection of the U number of different signals after radio 

frequency intermediate frequency downconversion able to cope with and exploit 
frequency selective multipath propagation, arbitrarily time-varying fading 
characteristics, non Line-Of-Sight communication. 

1 2. A method of communicating information from one location to another though 
1 5 a wavefield space using a multi-dimensional communication protocol comprising the 

steps of: 

receiving a message at an input terminal; 

encoding a plurality of signals associated with the message using an encoder, 
transmitting the plurality of signals associated with the message using a first 

20 plurality of antennas; 

receiving the plurality of signals associated with the message using a second 

plurality of antennas; and 

decoding the plurality of signals associated with the message to reconstruct the 
original information using a receiver, wherein the decoding step uses no advance 
25 knowledge of at least one of the dimensions of the multi-dimensional protocol. 

13. The method of claim 12, wherein: 

the decoding step includes the step of executing a statistical procedure to 
identify and decode the plurality of signals. 
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14. The method of claim 13, wherein: 

the decoding step includes the step of executing a generalized likelihood 
detection procedure to identify the plurality of signals and then executing a maximum 
likelihood detection procedure to decode the plurality of signals. 

5 15. The method of claim 14, wherein 

the decoding step includes the step of estimating the wavefield space using a 
virtual wavefield space that approximates the wavefield space. 

16. An apparatus for communicating information from one location to another 
though a wavefield space using a multi-dimensional communication protocol 

10 comprising: 

means for receiving a message at an input terminal; 
means for encoding a plurality of signals associated with the message; 
means for transmitting the plurality of signals associated with the message; 
means for receiving the plurality of signals associated with the message; and 
15 means for decoding the plurality of signals associated with the message to 

reconstruct the original information, wherein the means for decoding has no advance 
knowledge of at least one of the dimensions of the multi-dimensional protocol. 
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